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Surface  structures  and  electronic  states  of  metal  (M  =  Cu,  Ag  and  Au)  adsorbed  on  Ce02(l  1 1 )  surface  have 
been  investigated  using  first-principles  density-functional  theory  calculations  with  the  on-site  Coulomb 
interaction  taken  into  account.  Results  indicate  that  the  adsorption  is  accompanied  by  an  electron  charge 
transfer  between  neutral  metal  atom  and  neighboring  Ce4+  cation.  Charge  redistribution  with  electron 
localized  on  the  various  Ce  cations  neighboring  to  adsorbed  atom  is  observed  and  can  be  attributed  to  local 
structural  distortion  effects.  For  the  same  adsorption  site,  the  structural  configuration  with  the  electron 
localized  on  the  surface  Ce  cation  away  from  the  adsorbed  metal  atom  M  is  energetically  favorable.  The 
strongest  adsorption  interaction  for  Cu  and  Ag  is  located  at  3-fold  subsurface  O-top  site,  whereas  2-fold 
surface  O-bridge  site  is  the  most  stable  for  Au  adsorption.  The  order  on  the  adsorption  energy  follows  the 
sequence  of  Cu  >  Ag  >  Au,  indicating  that  Cu  atom  has  a  stronger  interaction  with  Ce02  (111)  surface  than 
Ag  and  Au.  Previous  controversies  about  the  localization  of  the  transferred  electron  and  the  valence  state 
of  the  adsorbed  Au  atom  have  been  successfully  clarified  via  analysis  of  the  local  structural  distortions. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Ceria-based  oxides  have  been  regarded  as  one  of  the  most 
important  components  in  many  catalytic  systems  due  to  their 
unique  properties  and  extensive  applications  as  active  supports  in 
metal-based  catalysts  for  oxidation  reactions  [1].  These  applica¬ 
tions  are  mainly  based  on  the  fact  that  ceria  can  act  as  an  oxygen 
reservoir  by  releasing/uptaking  oxygen  through  the  redox  pro¬ 
cesses  of  the  Ce4+/Ce3+  couple  [2-4]  and  the  strong  interactions 
between  metal  atom  and  ceria  support.  For  example,  increasing 
the  number  and  mobility  of  oxygen  vacancies  enhances  the  OSC  and 
correspondingly  improves  the  catalytic  activity  [5,6].  The  reducibil- 
ity  and  catalytic  activity  of  Ce02  are  significantly  enhanced  by  the 
presence  of  a  small  amount  of  transition  metal  atoms,  i.e.,  Cu,  Ag, 
and  Au.  However,  to  date,  basic  issues,  such  as  the  stable  adsorption 
sites  and  the  electronic  states  at  M/Ce02  (M  =  Cu,  Ag,  Au)  interfaces 
remain  open.  Especially  a  clarification  of  the  charge  state  of  Au  atom 
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on  stoichiometric  Ce02(l  1 1)  surface  is  urgently  needed,  which  is 
one  of  the  aims  of  the  present  work. 

Associated  with  this,  the  Cu/Ce02  system  has  been  identified  as 
an  excellent  catalyst  for  a  variety  of  important  reactions,  including 
direct  oxidation  of  hydrocarbon  fuels  in  solid  oxide  fuel  cells  (SOFC) 
[7-9],  NO  reduction  [10,11]  and  CO  oxidation  [12,13].  Previously, 
several  studies  indicated  that  the  Cu-O-Ce  structures  and  the  par¬ 
ticular  electronic  states  of  Cu  and  Ce  play  an  important  role  in  the 
Cu/Ce02  based  catalysts  [14,1 5].  The  presence  of  both  Cu2+  and  Cu1+ 
species  in  the  catalysts  has  been  found  by  several  research  groups 
[14-18].  In  a  recent  work,  Ye  et  al.  [19]  suggested  that  Cu  loading 
on  Ce02  not  only  provides  more  conductive  paths  for  electrons, 
but  also  connects  more  effective  reaction  sites,  causing  an  enhanc¬ 
ing  catalytic  activity  of  ceria.  Using  the  first-principles  calculations, 
Yang  et  al.  [20]  and  Branda  et  al.  [21  ]  investigated  the  properties  of 
Cu  adsorbed  on  the  Ce02(l  1 1)  surface,  and  their  results  demon¬ 
strate  the  existence  of  Cu1+  on  the  surface  of  Ce02,  although  the 
extra  electron  is  found  to  delocalize  on  three  surface  Ce  ions  which 
contradicts  the  common  sense. 

The  Ag/Ce02  system  is  less  studied  compared  with  the  cases 
of  Cu/Ce02  and  Au/Ce02.  Scire  et  al.  [22]  prepared  Ag/Ce02  sam¬ 
ple  and  assumed  that  its  high  catalytic  activity  in  the  methanol 
combustion  is  mainly  due  to  the  formation  of  smaller  Ag  nanopar¬ 
ticles.  The  same  conclusions  have  been  reported  by  Beier  et  al.  [23]. 
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They  suggested  that,  the  Ag  exhibits  a  similar  or  even  more  supe¬ 
rior  catalytic  effect  in  the  presence  of  Ce02,  compared  with  Pd 
and  Au  catalysts.  To  examine  the  reduction  processes  of  oxygen 
molecules  on  silver  cathode  surfaces,  Wang  et  al.  [24]  investigated 
the  monolayer  Ag  adsorbed  on  the  Ce02  surfaces,  implying  that 
oxygen  molecules  prefer  being  reduced  at  triple-phase  boundary 
(TPB). 

Au/Ce02  system  has  also  been  reported  to  be  a  prominent  cata¬ 
lyst  for  the  water-gas  shift  (WGS)  reaction  used  for  the  production 
of  hydrogen  in  the  chemical  industry  as  well  as  in  fuel  cell  appli¬ 
cations  [25-28].  The  origin  of  the  excellent  performance  of  the 
Au/Ce02  catalysts  has  been  extensively  studied  from  both  experi¬ 
mental  and  theoretical  points  of  view.  However,  the  disagreements 
appear  concerning  the  electronic  state  of  Au  adsorbed  on  the  sur¬ 
face  of  Ce02.  Fu  et  al.  [29,30]  and  Skoda  et  al.  [31]  suggested 
that  the  high  catalytic  activity  of  Au/Ce02  arises  from  positively 
charged  Au.  Using  the  first-principles  calculations,  Liu  et  al.  [32] 
also  suggested  that  the  high  catalytic  activity  is  mostly  caused  by 
Au  cations  instead  of  the  metallic  Au  particles.  A  different  conclu¬ 
sion  was  drawn  from  the  X-ray  absorption  fine-structure  (XAFS) 
spectroscopy  by  Wang  et  al.  [33,34],  indicating  that  Au  nanoparti¬ 
cles  remain  neutral.  Similarly,  Castellani  et  al.  [35]  investigated  the 
atomic  and  electronic  structures  of  the  Au  atom  on  ceria,  and  found 
that  Au  for  most  of  the  active  sites  remains  essentially  neutral. 
Zhang  et  al.  [36]  used  density  functional  theory  to  investigate  the 
nucleation  of  small  Au  clusters  (up  to  Aun ),  and  found  the  contact 
layer  Au  atoms  bound  with  surface  0  atoms  are  positively  charged, 
which  could  have  significant  implications  in  catalysis.  Recently, 
Branda  et  al.  [21]  suggested  that  the  description  of  Au  adsorption 
is  more  method  dependent  due  to  the  nearly  degeneracy  between 
the  cationic  and  neutral  solutions. 

Although  some  progresses  have  been  achieved,  there  are  still 
many  uncertainties  for  the  Cu/Ce02,  Ag/Ce02  as  well  as  Au/Ce02, 
such  as  the  valence  states  of  the  adsorbed  M  (M  =  Cu,  Ag,  and  Au), 
the  energetically  favorable  adsorption  sites  and  the  order  on  the 
adsorption  energy.  Especially,  it  is  a  scientifically  and  technolog¬ 
ically  important  issue  to  characterize  and  understand  the  charge 
transfer  and  the  subsequent  localization  of  electron  on  surface 
Ce  ions  upon  the  adsorption.  In  the  present  paper,  we  make  an 
attempt  to  find  out  the  inherent  adsorption  mechanism  of  Cu,  Ag, 
and  Au  on  Ce02(l  1 1 )  surface  using  the  first-principles  calculations 
based  on  the  density  functional  theory  (DFT),  following  our  previ¬ 
ous  studies  on  the  effects  of  M  (M  =  Mn,  Pr,  Sn,  Zr)  doping  on  the 
redox  thermodynamics  of  bulk  Ce02  [37].  We  have  comprehen¬ 
sively  investigated  the  surface  structures  and  electronic  states  of 
the  adsorbed  Cu,  Ag,  and  Au  on  Ce02(l  1 1 )  surface.  In  particular,  to 
the  best  of  our  knowledge,  the  effect  of  charge  distribution  upon  the 
adsorption  due  to  the  local  structural  distortions  is  firstly  proposed 
to  understand  the  interaction  between  the  supported  metal  atoms 
and  Ce02(l  1 1 )  surface.  We  also  discuss  the  present  disagreements 
on  the  valence  state  of  the  adsorbed  Au  atom.  Moreover,  the  effects 
of  the  exchange-correlation  functionals  and  lattice  parameters  on 
the  electronic  states  of  Au/Ce02(  1 1 1 )  interface  are  clearly  commu¬ 
nicated  in  this  paper. 


2.  Computational  details 

The  calculations  are  performed  using  the  Vienna  ab  initio  simu¬ 
lation  package  (VASP)  in  the  framework  of  the  projector  augmented 
wave  method  [38-40].  The  Perdew-Burke-Ernzerhof  (PBE)  func¬ 
tional  is  used  for  the  exchange  correlation  [41  ]  unless  otherwise 
specified.  Furthermore,  in  order  to  account  for  the  strong  on-site 
Coulomb  repulsion  among  the  Ce  4f  electrons,  a  Hubbard  param¬ 
eter  U  is  added  to  the  PBE  functional  employing  the  rotationally 
invariant  approach  formulated  by  Dudarev  et  al.  [42],  in  which  only 
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Fig.  1.  Possible  adsorption  sites  for  M  ( M  =  Cu ,  Ag,  Au)  on  the  (2  x  2)  Ce02(l  1 1) 
surface:  (a)  side  and  (b)  top  views.  Red  and  blue  spheres  represent  0  and  Ce  atoms, 
respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  the  article.) 

the  difference  (Ueff  =  U-J)  between  the  Coulomb  U  and  exchange  J 
parameters  enters.  Here,  we  choose  a  Ueff  value  (hereafter  referred 
to  simply  as  U)  with  the  value  of  4.5  eV  as  suggested  by  Fabris 
et  al.  [43].  The  value  removes  the  self-interaction  error  and  thus 
improves  the  description  of  the  correlation  effects.  In  particular,  the 
density  functional  theory  (DFT)  method  has  been  shown  suitable 
to  describe  the  Ce02  [44-47]. 

In  our  calculations  the  configurations  of  Ce  5s25p66s25d14f1,  O 
2s22p4,  Cu  3d104s1 ,  Ag  4d105s1 ,  and  Au  Sd^s1  are  treated  as  the 
valence  electrons.  We  use  a  plane-wave  cutoff  energy  of 400  eV  and 
2x2x1  Monkhorst-Pack  k-point  mesh  [48]  with  Gaussian  smear¬ 
ing  of  0.20  eV  [49]  for  the  structural  optimization  and  total  energy 
calculations.  Total  energies  are  converged  to  10-7  eV  unit  cell-1. 
The  Hellmann-Feynman  forces  are  converged  to  be  less  than 
0.005  eV  A-1.  For  a  more  accurate  calculation  of  the  electronic 
density  of  states  (DOS),  we  apply  4x4x1  k-point  mesh  and  the 
modified  tetrahedron  method  [50]. 

The  stoichiometric  (111)  surface  is  modeled  by  a  p(  2  x  2)  unit 
cell  with  twelve  atomic  layers  and  a  vacuum  layer  of  15  A  is  intro¬ 
duced  to  separate  the  slabs  along  the  z  direction.  The  x  and  y 
dimensions  of  the  supercells  are  fixed  at  the  calculated  bulk  values, 
but  the  ions  herein  are  allowed  to  relax.  The  present  PBE  +  U  calcula¬ 
tions  give  an  equilibrium  lattice  parameter  of  5.49  A,  which  is  close 
to  the  experimental  value  of  5.41  A  [51  ].  To  model  the  M/Ce02  sys¬ 
tems,  the  M  atom  is  put  at  various  active  sites  on  one  side  of  Ce02 
p( 2  x  2)  surfaces  as  depicted  below.  The  bottom  six  layers  are  fixed 
at  their  bulk  positions  to  mimic  the  bulk  and  remaining  layers  are 
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(a)  Od  (3-fold) 


(b)Ou  (O-top) 


(c)  Ou-Ou  (O-bridge)  (d)  Ce-Ce  (Ce-bridge) 


Fig.  2.  Structural  models  of  M  (M  =  Cu,  Ag,  Au)  adsorbed  on  (a)  Od  (3-fold),  (b)  Ou  (O-top),  (c)  Ou-Ou  (O-bridge)  and  (d)  Ce-Ce  (Ce-bridge)  sites.  Red,  blue  and  purple  spheres 
represent  0,  Ce  and  M  atoms,  respectively.  The  notations  of  INN  to  3NN  indicate  the  Ce  atoms  which  are  the  first,  second  and  third  nearest  neighbors  to  the  adsorbed  M 
atom,  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


allowed  to  relax.  Prior  to  VASP  optimization,  small  structural  dis¬ 
tortions  are  introduced  into  the  as-built  models  to  achieve  various 
localized  solutions.  These  require  slight  adjustments  of  distances 
between  a  Ce4+  ion  cores  and  nearest-neighboring  O  ion  cores  to 
be  consistent  with  those  of  a  Ce3+  ion.  The  adsorption  energy  of 
atom  M  on  the  Ce02  surface  is  calculated  by 

£ads  =  £ce02  +  %  “  EM/Ce02  > 

where  £Ce o2»  %  and  EM/Ce q2  correspond  to  the  calculated  total 
energies  of  the  clean  Ce02  surface,  adsorbed  metal  atom  M  and 
combined  M/Ce02  system,  respectively. 

3.  Results  and  discussion 

To  explore  the  various  possible  adsorption  models  of  M  (M= Cu, 
Ag,  Au)  on  the  Ce02(  1 1 1)  surface,  the  single  M(M  =  Cu,Ag,  Au)  atom 
is  placed  on  the  different  positions  derived  from  the  fold  struc¬ 
ture  of  the  surface  O  atom  and  subsurface  Ce  atom  as  shown  in 
Fig.  1,  including  a  three-fold  subsurface  O-top  (Od)  site,  a  surface 
O-top  (Ou)  site,  a  surface  O-bridge  (Ou-Ou)  site,  and  a  subsurface 
Ce-bridge  (Ce-Ce). 

Associated  with  the  adsorption  of  M  (M=Cu,  Ag,  Au)  atom,  a 
charge  transfer  between  M  and  Ce02(  1 1 1 )  surface  has  been  clearly 
observed,  causing  the  appearance  of  well  characterized  Ce3+  and 
M+  cations.  It  should  be  noted  that,  one  can  find  various  configu¬ 
rations  for  4f  electrons  to  locate  at  different  neighboring  Ce  ions 
to  the  M  (M  =  Cu,  Ag,  Au)  atom,  even  for  the  same  adsorption  site. 
As  an  example,  for  the  three-fold  Od  site  as  shown  in  Fig.  2a,  in 
the  outmost  cerium  plane,  one  can  find  three  nearest  equivalent  Ce 
cations  (INN)  and  one  second  nearest  Ce  cations  (2NN)  around  the 
adsorbed  metal  atom.  We  find  that  the  electron  can  actually  localize 
on  the  outmost  surface  Ce  ion  at  the  INN  or  2NN  to  the  metal  atom 
with  the  similar  adsorption  energies.  The  same  case  applies  for  the 
one-fold  Ou  site  and  Ce-Ce  bridge  site,  i.e.,  there  are  also  two  kinds 
of  charge  transfer  configurations  [see  Fig.  2b  and  d].  However,  for 
the  two-fold  Ou-Ou  site,  there  are  one  1 NN  Ce  ion,  two  2NN  Ce  ions 


and  one  3NN  Ce  ion  around  the  metal  atom,  therefore,  the  electron 
can  localize  on  the  INN,  2NN  or  3NN  Ce  ion  [see  Fig.  2c]. 

3.1.  Cu/Ce02(  1  1  1J 

Based  on  the  configurations  proposed  above,  we  calculate  iso¬ 
surfaces  of  excess  spin  density  for  Cu  adsorbed  on  Ce02(l  1 1) 
surface  as  shown  in  Fig.  3.  It  is  interesting  to  see  that,  depending 
on  the  local  atomic  configurations,  the  electron  transferred  from 
the  Cu  atom  can  localize  on  the  INN,  2NN  or  3NN  Ce  ion  to  the 
Cu.  Table  1  summarizes  the  adsorption  energies,  geometries  and 
electronic  properties  of  all  the  possible  adsorption  modes  of  Cu 
adsorbed  on  Ce02(l  1 1)  surface.  It  can  be  clearly  seen  that  all  the 
considered  adsorptions  involve  the  reduction  of  one  of  the  surface 
cerium  atoms,  as  indicated  by  magnetic  moments  between  0.82  and 
0.96 /xb,  which  are  close  to  the  ideal  localized  magnetic  moment 
of  1  /xB  per  Ce3+.  Furthermore,  the  local  environments  of  the  Ce 
atoms  differ  from  each  other  due  to  the  larger  size  of  the  Ce3+  ion 
compared  with  the  Ce4+  ion.  The  average  bond  lengths  between 
the  localized  Ce3+  ion  and  the  neighboring  O  ions  for  the  respec¬ 
tive  configurations  are  in  the  range  of  2.45-2.54  A,  which  are  larger 
than  the  theoretical  average  bond  length  (2.38  A)  of  Ce4+-0,  but 
are  comparable  to  the  experimental  Ce3+-0  bond  length  of  2.50  A 
[52].  Also,  this  is  in  reasonable  agreement  with  the  average  bond 
length  (2.47  A)  of  Ce3+-0  on  the  O-deficient  Ce02(l  1 1)  surface 
[52].  Actually,  the  average  bond  lengths  of  Ce4+-0  for  various  con¬ 
figurations  are  identical  to  each  others.  The  bond  lengths  between 
Cu  and  surface  O  atoms  are  in  the  range  of  1.73-2.07  A  with  the 
value  of  2.03  A  for  the  most  stable  configuration,  which  is  consistent 
with  the  experimental  value  of  1.86  A  in  cubic  Cu20  [53].  It  should 
be  mentioned,  when  the  adsorption  occurs  at  the  Ou-Ou  site,  the 
adsorption  energies  obey  the  sequence:  3NN  (conf.  1)>2NN  (conf. 
2)>1NN  (conf.  3).  For  the  Od  and  Ou  sites,  the  adsorption  at  the 
2NN  (conf.  1)  site  is  also  more  energetically  favorable  than  that  at 
the  INN  (conf.  2)  site.  Therefore,  we  may  assume  that  for  the  same 
adsorption  site,  the  excess  electrons  localize  preferentially  at  the  Ce 
ion  relatively  far  way  from  the  adsorbed  Cu.  Note  that,  the  systems 
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Od  (conf.  2,  INN)  Ou  (conf.  1,  2NN) 


[1 


[0  1  -1] 

Od  (conf.  1,  2NN) 


[0  1  -1] 

Ou  (conf.  2,  INN)  0„-0„  (conf  1, 3NN)  Ou-Ou  (conf.  2,  2NN) 


[0  1  -1] 

Ou-Ou  (conf.  3,  INN)  Ce-Ce(lNN) 


Fig.  3.  Isosurfaces  of  excess  spin  charge  density  for  Cu  adsorbed  on  different  sites  of  Ce02(l  1 1 )  surface.  The  labeling  corresponds  to  the  description  of  configuration  and  the 
location  of  the  Ce3+  ion  in  Table  1  and  Fig.  2. 


have  to  be  distorted  slightly  beforehand  to  obtain  these  configura¬ 
tions,  causing  the  atomic  environment  of  the  desired  initial  position 
Ce4+  ions,  to  a  large  extent,  to  be  consistent  with  that  of  the  Ce3+. 
To  our  knowledge,  the  phenomenon  that  the  slight  structural  dis¬ 
tortion  can  cause  the  complete  oxidization  of  adsorbed  metal  atom 
M  in  M/Ce02  interfaces  has  never  been  reported  in  the  literatures, 
although  similar  one  has  been  observed  from  the  investigations  of 
the  O-vacancy  on  the  Ce02(l  1 1)  surface  by  Ganduglia-Pirovano 
et  al.  [52]  and  Li  et  al.  [54],  which  indicated  that  there  exist 
multiple  configurations  of  the  two  excess  electrons  localization 
associated  with  the  formation  of  an  O-vacancy  on  the  Ce02(l  1 1) 
surface. 

From  our  calculations,  the  Cu  adsorption  at  the  three-fold  Od  site 
is  the  most  energetically  favorable  with  the  adsorption  energy  of 
2.62  eV,  which  is  consistent  with  the  results  of  Branda  et  al.  [21  ]  and 
Yang  et  al.  [20].  However,  it  should  be  emphasized  that  all  calcula¬ 
tions  on  the  Cu  adsorption  at  the  Od,  Ou-Ou  and  Ou  sites  by  Branda 
et  al.  [21  ]  and  Yang  et  al.  [20]  indicate  that  electron  released  from 
neutral  metal  atom  delocalizes  on  all  three  nearest  neighboring 
(NN)  Ce  ions,  which  is  different  from  our  conclusion,  i.e.,  the  trans¬ 
ferred  charge  should  localize  on  one  2NN  Ce  ion  with  the  magnetic 
moment  of  0.96  /xB.  In  fact,  the  electron  transfer  from  Cu  to  a  single 
Ce  ion  can  be  applied  for  all  the  configurations  in  our  work.  We 
also  find  the  delocalized  structures  obtained  are  less  stable  than 
the  localized  ones  and  the  corresponding  Ce4+  is  not  completely 
reduced  to  Ce3+.  Therefore,  we  argue  that  the  localized  solution  is 
more  reasonable.  A  similar  conclusion  has  also  been  drawn  from 
our  previous  investigation  on  the  O-vacancy  in  bulk  Ce02,  in  which 
the  two  excess  electrons  upon  the  formation  of  the  O-vacancy  pre¬ 
fer  localizing  on  two  NN  Ce  ions  instead  of  delocalizing  on  four  NN 
Ce  ions  [37]. 


Since  the  adsorption  is  accompanied  by  a  charge  transfer  from 
Cu  to  an  unoccupied  Ce  4f  level,  it  is  significant  to  understand 
the  electronic  properties  of  the  adatom  Cu  and  substrate  atoms. 
Here  we  take  the  most  energetically  favorable  Od  (conf.l)  model 
for  an  example  and  analyze  the  respective  partial  density  of  states 
(PDOS)  of  Cu  atom,  pure  Ce02(l  11),  and  Cu/Ce02(l  1 1)  as  shown 
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Fig.  4.  Partial  density  of  states  (PDOS)  of  Cu  atom,  pure  Ce02(l  1 1)  and  Cu  adsorbed 
on  Od  site  (conf-1,  2NN).  The  zero  point  of  the  energy  axis  corresponds  to  the  Fermi 
level,  whereas  the  upper  and  lower  panels  represent  the  spin  up  and  spin  down 
channels,  respectively. 


32 


Y.  Tang  et  al.  /  Journal  of  Power  Sources  197  ( 2012)28-37 


Table  1 

Adsorption  energies,  geometries  and  electronic  properties  for  the  adsorption  configurations  of  Cu/Ce02(l  11).  d(ai-o)  and  d(cu-ce)  are  distances  from  the  adsorbed  Cu  to 
surface  0  and  Ce  atoms,  respectively.  /xCe  and  /x cu  are  the  magnetic  moments  of  reduced  Ce  (in  bold)  and  oxidized  Cu,  respectively.  d(Ce-o)  corresponds  to  the  average  Ce3+-0 
bond  length.  Previous  results  are  also  listed  for  comparison. 


Site 

fads  (eV) 

d(c u-o)  (A) 

d(Cu-Ce)  (A) 

Mce  (Mb) 

Mcu  (Mb) 

d(Ce-O)  (A) 

Reference 

3x2.88 

0.00 

Od  (conf.l,  2NN) 

2.62 

3  x  2.03 

1  x  4.75 

0.96 

0.00 

2.54 

Present  work 

1  x  2.04 

2x2.87 

0.01 

2.10 

2  x  2.07 

1  x  2.89 

0.86 

0.00 

2.46 

Present  work 

Od  (conf.2,  INN) 

1  x  4.83 

0.10 

1.86 

1  x  2.01 

2x2.78 

0.16,  0.19 

0.00 

[21] 

2  x  2.04 

1  x  2.79 

0.34 

VASP-(L5)L5//G3a 

2.68 

1  x  2.03 

1  x  2.88 

2x3.70 

0.00 

[20] 

Ou  (conf.l,  2NN) 

1.71 

1  x  1.77 

1  x  3.81 

0.00 

0.01 

2.48 

Present  work 

1  x  5.22 

0.89 

2  x  3.60 

0.00 

1.63 

1  x  1.78 

1  x  3.88 

0.91 

0.02 

2.46 

Present  work 

1  x  4.95 

0.00 

Ou  (conf.2,  INN) 

1  x  3.57 

0.18 

[21] 

1.50 

1  x  1.73 

1  x  3.61 

0.16 

0.08 

VASP-(L5)L5//G3a 

1  x  3.65 

0.22 

1.63 

[20] 

Ou-Ou  (conf.l, 

3NN) 

2.45 

2x1.94 

1  x  2.77 

2  x  3.02 

1  x  4.66 

0.00 

0.00 

0.96 

0.00 

2.52 

Present  work 

1  x  2.70 

0.00 

Ou-Ou  (conf.2, 

2.17 

1  x  1.90 

1  x  3.19 

0.94 

0.00 

2.48 

Present  work 

2NN) 

1  x  1.93 

1  x  3.30 

0.00 

1  x  4.51 

0.00 

1  x  2.77 

0.90 

Ou-Ou  (conf.3, 

2.06 

2x1.93 

2  x  3.07 

1  x  4.79 

0.00 

0.00 

0.00 

2.51 

Present  work 

INN) 

1.61 

2  x  1.89 

1  x  2.61 

0.28 

0.00 

[21] 

2x3.13 

0.30,  0.45 

VASP-(L5)L5//G3a 

2x3.35 

0.02 

Ce-Ce 

1.44 

1  x  1.80 

1  x  4.17 

0.86 

0.03 

2.45 

Present  work 

1  x  5.72 

0.00 

a  The  notation  of  (L5)L5//G3  indicates  that  lattice  parameter  and  geometry  optimization  are  carried  out  at  the  LDA  +  U  (U=  5  eV)  level  whereas  the  energy  and  electronic 
properties  arise  from  single  point  calculations  at  the  GGA  +  U  (U= 3  eV)  level. 


in  Fig.  4.  It  can  be  seen  from  Fig.  4a  that  for  the  Cu  atom,  there 
is  an  asymmetric  occupied  peak  corresponding  to  the  Cu  4s  state. 
However,  the  adsorption  of  Cu  at  the  Ce02(l  1 1 )  surface  obviously 
gives  rise  to  the  disappearance  of  the  occupied  Cu  4s  gap  state  as 
shown  in  Fig.  4b,  indicating  that  the  single  Cu  4s  electron  has  been 
released  in  Cu/Ce02  system.  Moreover,  there  is  a  new  occupied  peak 
corresponding  to  the  Ce4f  states  for  the  Cu/Ce02(l  1 1),  further  con¬ 
firming  that  the  electron  indeed  localizes  on  the  2NN  Ce  ion,  thus 
making  Ce4+  and  Cu  reduced  and  oxidized  to  Ce3+  and  Cu+,  respec¬ 
tively.  This  is  in  good  agreement  with  the  previous  experimental 
results,  which  indicate  that  the  adsorbed  Cu  on  Ce02(l  1 1 )  surface 
would  be  oxidized  to  Cu+  [14,55,56],  which  is  also  consistent  with 
suggestion  of  Hornes  et  al.  [57]  that  a  high  catalytic  activity  of  CO 
oxidation  of  Ce02/CuO  system  are  related  to  interfacial  Cu+  species 
generated  during  the  reaction. 

3.2.  Ag/Ce02( 1  1  1) 

Following  the  same  methodology,  we  extend  our  studies  to 
the  Ag/Ce02(l  1 1)  system,  as  shown  in  Table  2.  We  also  find  the 
adsorption  on  the  three-fold  Od  site  is  the  most  energetically  favor¬ 
able  and  correspondingly,  the  single  electron  transfers  from  the 
adsorbed  Ag  neutral  atom  to  the  2NN  Ce4+  cation  although  there 
are  some  variances  from  Cu/Ce02(l  11).  Note  that,  both  Ou-Ou  and 
Ce-Ce  models  in  Ag/Ce02  system  are  unstable  and  converged  to  the 
Od  one  upon  the  structural  optimization.  This  is  in  a  good  agree¬ 
ment  with  the  results  of  Branda  et  al.  [21].  Therefore,  it  can  be 
concluded  that  only  Od  and  Ou  adsorption  configurations  exist  for 


the  Ag/Ce02(l  1 1)  system  by  using  the  density  functional  theory 
study. 

It  is  also  found  from  Table  2  that  the  transferred  electron  for  all 
the  calculated  models  localizes  on  a  single  Ce  ion.  Similar  to  the  case 
of  Cu/Ce02(l  1 1)  system,  there  exist  the  different  Ce  4f  electron 
localization  configurations  for  the  Od  and  Ou  adsorption  sites,  and 
the  electron  preferentially  localizes  on  the  surface  Ce  ion  which  is 
relatively  far  away  from  the  adsorbed  Ag,  which  from  results  in  Ref. 
[21] 

Magnetic  moments  of  the  localized  Ce  ions  for  all  the  configura¬ 
tions  are  in  the  range  of  0.80-0.96  /xB,  which  are  close  to  1  /xB/Ce3+. 
Moreover,  the  calculated  average  bond  lengths  of  Ag+-0  in  respec¬ 
tive  configurations  are  in  the  range  of  2.01  -2.47  A  with  three  Ag+-0 
bond  lengths  of  2.36, 2.36  and  2.38  A  for  the  most  stable  configura¬ 
tion,  which  are  slightly  larger  than  the  experimental  value  of  2.05  A 
in  cubic  Ag20  [58].  This  can  be  understood  from  the  fact  that  three 
and  two-fold  Ag-0  bonds  occur  in  the  Ag/Ce02(l  1 1)  (Od  adsorp¬ 
tion  configuration)  and  cubic  Ag20,  respectively.  These  analyses 
demonstrate  that  Ag  atom  is  indeed  oxidized  to  Ag+.  This  is  con¬ 
sistent  with  the  result  presented  by  Branda  et  al.  [21],  although 
they  did  not  observe  the  multiple  configurations  with  the  excess 
electron  localized  at  different  neighboring  Ce  ions.  The  similar  con¬ 
clusion  has  been  drawn  by  Wang  et  al.  [24],  who  investigated  a 
Ce02(l  1 1)  supported  atomic  layer  of  silver  using  a  (2x2)  slab 
model  with  four  Ag  atoms.  Results  indicated  that  there  is  0.1 4  elec¬ 
tron  on  the  adsorbed  monolayer  Ag,  which  reveals  that  all  four  Ag 
atoms  are  nearly  oxidized  to  Ag+.  Experimentally,  Farmer  et  al.  [59] 
studied  the  Ag  adsorption  on  reduced  Ce02(l  1 1)  thin  films  using 
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Table  2 

Adsorption  energies,  geometries  and  electronic  properties  for  the  adsorption  configurations  of  Ag/Ce02(l  11).  d(Ag_ o)  and  d(Ag_ ce)  are  distances  from  the  adsorbed  Ag  to 
surface  0  and  Ce  atoms,  respectively.  /zCe  and  /xAg  are  the  magnetic  moments  of  reduced  Ce  (in  bold)  and  oxidized  Ag,  respectively.  d(Ce- o)  corresponds  to  the  average  Ce3+-0 
bond  length.  Previous  results  are  also  listed  for  comparison. 


Site 

fads  (eV) 

d(Ag-o)  (A) 

d(Ag-Ce)  (A) 

MCe  (P-b) 

^Ag  (P-b) 

d(Ce- 0)  (A) 

Reference 

2  x  2.36 

1x3.18 

0.01 

Od  (conf.2,  2NN) 

1.42 

1  x  2.38 

2x3.19 

0.00 

0.00 

2.50 

Present  work 

1  x  4.94 

0.96 

2  x  2.47 

1  x  3.17 

0.90 

Od  (conf.l,  INN) 

1.11 

1  x  2.34 

2x3.15 

1  x  4.98 

0.00 

0.01 

0.00 

2.47 

Present  work 

1.00 

3  x  2.28 

3x3.07 

0.89 

0.00 

[21] 

2x3.86 

0.01 

Ou  (conf.l,  2NN) 

0.86 

1  x  2.07 

1  x  4.00 

0.00 

0.02 

2.47 

Present  work 

1  x  5.31 

0.82 

2x3.76 

0.00 

0.83 

1  x  2.08 

1  x  4.11 

0.80 

0.03 

2.45 

Present  work 

Ou  (conf.2,  INN) 

1  x  5.02 

1  x  3.65 

0.00 

0.56 

1  x  2.01 

1  x  3.75 

1  x  3.81 

0.64 

0.07 

[21] 

adsorption  microcalorimetry  and  found  that  the  heat  of  adsorption 
of  Ag  vapor  is  2.07  eV  (200  kj  mol-1 ),  which  is  much  larger  than  the 
present  value  of  1.42  eV.  This  may  be  attributed  to  the  existence  of 
O-vacancies  on  reduced  Ce02(l  11). 

3.3.  Au/Ce02(l  11) 

In  spite  of  much  attention  paid  on  Au/Ce02  system,  both  the 
oxidizations  states  of  the  adatom  Au  and  energetically  favor¬ 
able  adsorption  configuration  are  still  pending  issues.  The  present 
PBE  +  U  calculations  indicate  that  Au  located  at  the  bridge-site  has 
the  largest  adsorption  energy  which  is  consistent  with  the  one  by 
Hernandez  et  al.  [60]  and  Camellone  et  al.  [61],  Zhang  et  al.  [62] 
and  Branda  et  al.  [21  ]  However,  Liu  et  al.  [32],  Castellani  et  al.  [35] 
and  Chen  et  al.  [63]  stated  that  Au  adsorption  at  the  O-top  of  the 
Ce02(l  1 1)  surface  is  the  most  stable.  The  disagreement  lies  in  the 
choice  of  the  ultrasoft  pseudopotentials,  the  value  of  the  U  or  cut¬ 
off  energy  used  to  describe  the  valence  electrons  [21,35].  Table  3 
summarizes  the  adsorption  energies,  relevant  bond  lengths  and 
magnetic  moments  for  a  single  Au  atom  on  the  Ce02(l  1 1 )  surface 
together  with  the  theoretical  and  experimental  results  reported  in 
the  literatures  [21,60-63].  The  most  stable  adsorption  configura¬ 
tion  is  the  Ou-Ou  model  with  the  electron  localized  on  the  3NN  Ce 
ion  to  the  Au.  The  adsorption  energy  and  Au-0  bond  length  in  the 
present  work  are  0.99  eV  and  2.15  A,  respectively,  which  are  very 
close  to  those  (1.15  eV  and  2.1 5  A)  of  Hernandez  et  al.  [60]  using 
the  GGA  +  U  formulation  with  the  U  value  of  5.0  eV. 

Accordingly,  it  can  be  seen  from  Table  3  that  there  are  also 
different  configurations  with  the  excess  electron  localized  on  the 
different  Ce  ions.  In  other  words,  similar  to  the  cases  of  the 
Cu/Ce02(lll)  and  Ag/Ce02(l  1 1)  systems,  the  excess  electron 
derived  from  the  Au  adsorption  at  the  Ou,  Od  and  Ou-Ou  sites  can 
condense  into  the  localized  4f-level  traps  of  the  INN,  2NN  or  3NN 
Ce  ion.  Moreover,  for  the  same  adsorption  site,  the  longer  the  dis¬ 
tance  between  the  surface  Ce  ion  with  localized  electron  and  the 
adsorbed  Au  atom  is,  the  more  energetically  favorable  the  config¬ 
uration  becomes. 

Moreover,  Zhang  et  al.  [62]  suggested  that  there  is  a  distribution 
of  Au°  (hollow  site)  or  Au+(top  and  bridge  site)  atomic  species  at  the 
regular  sites  of  the  Ce02(l  1 1 )  surface  and  Baron  et  al.  [64]  believed 
that  there  exists  partially  charged  Au5+  on  Ce02(l  11).  Recently, 
Branda  et  al.  [21,65]  argued  that  the  description  of  Au/Ce02(l  1 1) 
is  more  method  dependent  and  there  exist  two  well-defined  states 
with  similar  energy  but  exhibit  different  electronic  structures  and 
local  geometries,  i.e.,  one  with  fully  oxidized  Au  and  the  other  with 
neutral  Au.  However,  our  calculations  indicate  that,  independent 


of  the  location  of  Au  at  the  surface  site,  the  adsorption  of  Au  atom 
on  Ce02(l  1 1)  surface  is  always  accompanied  by  the  one-electron 
charge  transfer  from  Au  atom  to  a  neighboring  Ce  ions,  which 
is  independent  on  the  adsorption  sites.  As  an  example  shown  in 
Table  3,  in  the  case  of  adsorption  on  the  Ou-Ou  site  with  the  elec¬ 
trons  transferred  to  the  3NN  Ce  ion,  the  magnetic  moment  of  Ce3+ 
is  0.97  ii b  and  the  Au  atom  is  oxidized  to  Au+  ion. 

To  further  analyze  the  interaction  between  the  adsorbed  Au 
atom  and  Ce02(l  1 1 )  surface,  we  examine  the  degree  of  the  charge 
transfer  between  Au  and  the  Ce02(l  1 1)  surface.  Bader  charge 
based  analysis  [66]  predicts  that  the  adsorbed  Au  loses  about  1.3e, 
revealing  a  completely  oxidization  process.  The  analysis  of  PDOS  for 
Au/Ce02(l  1 1)  provides  a  further  confirmation.  It  can  see  seen  from 
Fig.  5  that  its  6s  state  upon  the  Au  atom  adsorption  becomes  unoc¬ 
cupied  and  a  new  peak  emerges  between  O  2p  valance  band  top 
and  unoccupied  Ce  4f  band  bottom.  In  other  words,  the  adsorbed 
Au  is  oxidized  whereas  the  3NN  Ce4+  is  reduced  to  Ce3+.  Actually, 
one  can  observe  always  a  single  neighboring  completely  reduced  Ce 
ion  for  all  the  considered  adsorption  configurations.  Experimental 
results  also  support  the  present  conclusion.  For  example,  Fu  et  al. 
[30]  found  Au+  by  using  the  X-ray  photoelectron  spectroscopy  (XPS) 
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Fig.  5.  Partial  density  of  states  (PDOS)  of  Au  atom,  pure  Ce02(l  1 1 )  and  Au  adsorbed 
on  Ou-Ou  site  (conf-1,  3NN).  The  zero  point  of  the  energy  axis  corresponds  to  the 
Fermi  level,  whereas  the  upper  and  lower  panels  represent  the  spin  up  and  spin 
down  channels,  respectively. 
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Table  3 

Adsorption  energies,  geometries  and  electronic  properties  for  the  adsorption  configurations  of  Au/Ce02(l  11).  d(Au-o)  and  d(Au-ce)  are  distances  from  the  adsorbed  Au  to 
surface  O  and  Ce  atoms,  respectively.  /zce  and  /xAu  are  the  magnetic  moments  of  reduced  Ce  (in  bold)  and  oxidized  Au,  respectively.  d(Ce- o)  corresponds  to  the  average  Ce3+-0 
bond  length.  Previous  results  are  also  listed  for  comparison. 


Site 

fads  (eV) 

d(Au-o)  (A) 

d(Au-Ce)  (A) 

lice  (p-b) 

Mau  (p-b) 

d(Ce-O)  (A) 

Reference 

2  x  2.35 

1  x  3.11 

0.01 

Od  (conf.l,  2NN) 

0.69 

1  x  2.36 

2x3.12 

0.00 

0.02 

2.49 

Present  work 

1  x  4.89 

0.91 

1  x  2.22 

2x3.15 

0.01 

0.67 

2  x  2.62 

1  x  3.19 

0.84 

0.02 

2.47 

Present  work 

1  x  4.98 

0.00 

Od  (conf.l,  INN) 

0.61 

3  x  2.76 

3x3.27 

1  x  5.06 

0.10 

0.02 

[60] 

1  x  2.56 

1  x  3.05 

0.51 

1  x  2.58 

1  x  3.08 

0.14 

0.25 

[21] 

1  x  2.63 

1  x  3.10 

0.66 

2x3.80 

0.01 

[62] 

0.90 

1  x  2.00 

1  x  3.95 

0.01 

0.08 

2.46 

Present  work 

1  x  5.26 

0.77 

2x3.85 

0.00 

[60] 

0.96 

1  x  1.96 

1  x  3.93 

0.00 

1  x  5.32 

0.96 

0.66 

1  x  2.07 

2  x  3.62 

1  x  3.63 

<0.1 

0.31 

[21] 

Ou  (conf.l,  2NN) 

1.26 

[32] 

1.04 

1  x  2.00 

[61] 

0.96 

[62] 

0.78 

1  x  2.18 

[63] 

SIESTA-PBE 

0.63 

1  x  2.18 

[63] 

SIESTA-BLYP 

0.88 

1  x  2.23 

1  x  3.69 

0.00 

[63] 

VASP 

0.85 

1  x  2.02 

1  x  3.70 

1  x  4.03 

0.00 

0.62 

0.15 

2.42 

Present  work 

Ou  (conf.2,  INN) 

1  x  5.01 

0.00 

2x3.50 

0.00 

0.92 

2.06 

1  x  4.22 

0.83 

[60] 

1  x  4.60 

0.00 

1  x  2.85 

0.00 

0.99 

2x2.15 

2x2.31 

0.00 

0.00 

2.49 

Present  work 

1  x  4.62 

0.97 

Ou-Ou  (conf.l, 
3NN) 

1.15 

2x2.15 

1  x  2.85 

2  x  2.32 

1  x  4.63 

0.00 

0.00 

0.97 

[60] 

1.18 

1  x  4.76 

[61] 

1.04 

1  x  2.14 

1  x  2.85 

2  x  3.35 

1  x  2.83 

0.00 

[62] 

0.86 

1  x  2.21 

2  x  3.38 

0.95 

0.00 

2.49 

Present  work 

1  x  4.67 

0.00 

Ou-Ou  (conf.2, 

1  x  2.15 

1  x  2.83 

0.00 

2NN) 

1.00 

1  x  2.21 

1  x  3.37 

0.96 

[60] 

1  x  3.40 

0.00 

1  x  4.66 

0.00 

1.25 

1  x  3.34 

1  x  2.87 

0.91 

0.70 

2  x  2.20 

2  x  3.29 

0.01 

0.00 

2.47 

Present  work 

1  x  4.76 

0.03 

Ou-Ou  (conf.3, 

1  x  2.87 

0.95 

INN) 

0.86 

2x2.21 

2  x  3.30 

0.00 

[60] 

1  x  4.75 

0.00 

0.43 

2x2.15 

1  x  2.76 

2x3.25 

2  x  0.47 

0.00 

[21] 

2  x  3.63 

0.00 

Ce-Ce 

0.85 

1  x  2.03 

1  x  4.13 

0.68 

0.00 

2.43 

Present  work 

1  x  5.67 

0.00 

technique.  They  also  suggested  the  catalytic  activity  was  affected 
by  the  ionic  state  Au  instead  of  metallic  Au. 

In  the  following  we  will  focus  on  the  effects  of  exchange- 
correlation  functional  and  lattice  parameters  on  the  oxidization 
states  of  Au.  Here,  following  the  calculation  procedures  in  the 
present  work,  we  examine  the  oxidization  states  of  Au  for  the  top 


Ou  site  model  using  both  GGA  +  H  {U=  5.0)  and  LDA  +  H  {U=  5.0) 
with  both  cell  parameters  a0  =  5.40  and  5.49  A  so  as  to  make  the 
overall  comparison  with  the  results  of  Branda  et  al.  [65],  as  shown 
in  Table  4.  Both  GGA  +  U  (U=  5.0)  and  LDA  +  U  (U=  5.0)  calculations 
for  the  case  of  a0  =  5.49  A  inevitably  lead  to  the  same  conclusions. 
That  is,  the  adatom  Au  is  almost  oxidized  and  one  surface  Ce4+ 
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Table  4 

Adsorption  energies,  geometries  and  electronic  properties  for  the  adsorption  configurations  of  Au/Ce02(l  1 1)  using  GGA  +  U  and  LDA  +  U  methods.  d(Au-o)  and  d(Au-ce)  are 
distances  from  the  adsorbed  Au  to  surface  0  and  Ce  atoms,  respectively.  /xCe  and  /xAu  are  the  magnetic  moments  of  reduced  Ce  (in  bold)  and  oxidized  Au,  respectively.  d(Ce- o) 
corresponds  to  the  average  Ce3+-0  bond  length.  Previous  results  are  also  listed  for  comparison. 


Site 

fads  (eV) 

d( Au-O)  (A) 

d(Au-Ce)  (A) 

Mce  (Mb) 

Mau  (Mb) 

d(Ce-O)  (A) 

Reference 

GGA  +  U(a0  =  5.49,1/=  5.0) 

1  X  3.81 

0.00 

1  x  3.82 

0.00 

Conf.  1  (2NN) 

1.09 

1  X  1.98 

1  x  3.95 

0.01 

0.02 

2.48 

Present  work 

1  x  5.27 

0.87 

1  x  3.78 

0.00 

Conf.  2  (INN) 

0.99 

1  x  1.99 

1  x  3.79 

1  x  3.95 

0.00 

0.82 

0.06 

2.45 

Present  work 

1  x  5.18 

0.00 

LDA  +  U  (a0  =  5.49,  U=  5.0) 

2  x  3.68 

0.00 

Conf.  1  (2NN) 

1.57 

1  x  1.94 

1  x  3.84 

0.00 

0.03 

2.44 

Present  work 

1  x  5.18 

0.78 

2x3.56 

0.00 

Conf.  2  (INN) 

1.51 

1  x  1.96 

1  x  3.93 

0.57 

0.09 

2.40 

Present  work 

1  x  4.90 

0.00 

GGA  +  U(a0  =  5.49,1/=  5.0) 

0.96 

0.94a 

0.07 

[65] 

LDA  +  U  (a0  =  5.49,  U=5.0) 
GGA +  U(a0  =  5.40,  U=5.0) 

1.37 

0.93a 

0.07 

[65] 

1  x  3.81 

0.01 

Conf.  1  (2NN) 

0.91 

1  x  1.99 

1  x  3.82 

1  x  3.93 

0.01 

0.01 

0.06 

2.45 

Present  work 

1  x  5.23 

0.81 

1  x  3.67 

0.00 

Conf.  2  (INN) 

0.80 

1  x  2.03 

1  x  3.68 

1  x  4.06 

0.00 

0.68 

0.12 

2.41 

Present  work 

1  x  4.92 

0.01 

LDA  +  U  (a0  =  5.40,  U=5.0) 

2  x  3.66 

0.01 

Conf.  1  (2NN) 

1.54 

1  x  1.96 

1  x  3.82 

0.01 

0.10 

2.40 

Present  work 

1  x  5.12 

0.70 

2  x  3.53 

0.01 

Conf.  2  (INN) 

1.50 

1  x  1.99 

1  x  3.91 

0.53 

0.19 

2.36 

Present  work 

1  x  4.86 

0.01 

GGA  +  U(a0  =  5.40,  U=  5.0) 

0.61 

0.06-0.10b 

0.33 

[65] 

LDA  +  U  (a0  =  5.40,  U=5.0) 

1.27 

0.07-0.10b 

0.32 

[65] 

a  Fully  localized  on  one  Ce  atom,  2NN  to  the  O-Au  bond. 
b  Delocalized  over  four  Ce  atoms. 
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Fig.  6.  Partial  density  of  states  (PDOS)  of  Au  atom,  pure  Ce02(l  1 1 )  and  Au  adsorbed 
on  Ou  site  (conf-1,  2NN).  The  calculation  are  performed  with  GGA  +  U  (U=  5.0), 
ao  =  5.40  A.  The  zero  point  of  the  energy  axis  corresponds  to  the  Fermi  level,  whereas 
the  upper  and  lower  panels  represent  the  spin  up  and  spin  down  channels,  respec¬ 
tively. 


ion  is  reduced,  which  is  in  good  agreement  with  the  conclusion 
in  Ref.  [65].  As  for  the  case  of  a0  =  5.40A,  the  magnetic  moments 
of  Ce  ion  using  both  GGA  +  U  (Ll=5.0)  and  LDA  +  U  (U=  5.0)  calcula¬ 
tions  are  found  to  be  0.53-0.81  /xB,  indicating  Au  is  also  oxidized. 
This  disagrees  with  the  suggestion  of  Branda  et  al.  [65],  in  which 
the  magnetic  moment  of  Ce  ion  around  the  adatom  Au  is  small 
(/xCe  <  0.1 /xB).  We  notice  that  Branda  et  al.  [65]  did  not  supply 
the  electronic  structure  analysis.  In  order  to  further  reveal  the  oxi¬ 
dized  state  of  the  adatom  Au  and  make  comparison  with  the  PBE  +  U 
calculations,  we  analyze  the  PDOS  for  GGA  +  U(U=  5.0),  a0  =  5.40  A 
with  electrons  localized  on  the  2NN  Ce  ions,  as  shown  in  Fig.  6.  We 
can  see  that  the  adsorption  of  Au  atom  causes  the  occurrence  of 
a  new  occupied  peak  corresponding  to  localized  Ce  4f  states  [see 
Fig.  6b].  Meanwhile,  the  occupied  Au  6s  states  in  Fig.  6a  are  shifted 
above  the  Fermi  level,  indicating  that  Au  indeed  loses  electron. 
This  is  in  good  agreement  with  the  result  of  PDOS  using  PBE  +  D 
method.  Therefore,  we  believe  that,  structural  distortion  instead  of 
the  reasonable  choice  of  exchange-correlation  functional  and  lat¬ 
tice  parameters  as  previously  suggested  [65],  is  more  important 
to  achieve  the  oxidization  of  Au  and  a  localized  solution.  Tests  are 
performed  for  Cu/Ce02(l  1 1)  and  Ag/Ce02(l  1 1)  systems  by  vary¬ 
ing  exchange-correlation  functional  and  lattice  parameters  and 
reveal  no  dependence  of  our  conclusions  on  the  choice  of  adsorp¬ 
tion  species.  In  addition,  both  Branda  et  al.  [65]  and  we  report  that 
adsorption  energies  using  LDA  +  U  calculations  are  larger  than  those 
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using  GGA  +  U  and  PBE  +  U  ones  which  is  largely  because  the  LDA 
somewhat  overestimates  surface  energies  [67]. 


4.  Discussion 

Based  on  the  above  results,  we  can  conclude  in  the  cases  of  M 
(M  =  Cu,  Ag,  Au)  adsorbed  on  Ce02(l  1 1)  surface,  Cu  and  Ag  pre¬ 
fer  to  be  adsorbed  on  the  Od  site,  whereas  Au  is  at  the  Ou-Ou 
site.  We  notice  that  Yang  et  al.  [68,69],  who  studied  Pd  and  Pt 
adsorbed  on  Ce02(l  11),  found  that  Pd  and  Pt  prefer  to  be  adsorbed 
on  the  O-bridge  (Ou-Ou)  and  3-fold  Od  sites,  respectively.  How¬ 
ever,  the  electron  is  found  to  be  delocalized  around  the  adsorption 
site  and  adatoms  Pd  and  Pt  are  only  partially  oxidized.  Recently 
Bruix  et  al.  [70]  predicted  that  the  most  favorable  interaction 
for  the  Pt/Ce02(l  1 1)  corresponds  to  adsorption  on  top  of  a  site 
bridging  two  nearest-neighbor  surface  O  atoms  in  a  closed-shell 
electronic  configuration  with  an  essentially  neutral  or  slightly  pos¬ 
itively  charged  Pt.  Nevertheless,  for  this  site  another  slightly  higher 
energy  solution  is  found  with  two  unpaired  electrons  in  the  unit 
cell  and  the  adsorption  represents  a  completely  different  physical 
picture,  where  the  adatom  Pt  has  been  oxidized  to  Pt+  with  the 
concomitant  reduction  of  a  single  Ce4+  to  Ce3+.  Besides,  Chafi  et  al. 
[71]  investigated  the  Ni-Ce02  interaction  using  the  GGA  scheme 
and  found  that  the  Ou  top  site  is  the  most  stable.  We  can  say  that 
the  electronic  states  of  Pt  (Pd,  Ni)/Ce02  systems  mentioned  above 
are  somewhat  contentious,  which  should  be  reasonably  character¬ 
ized  via  analysis  of  the  local  structural  distortions.  In  other  words, 
for  the  purpose  of  achieving  the  energetically  stable  structures, 
the  interaction  of  metal  adatom  with  stoichiometric  Ce02(l  1 1) 
surface  will  result  in  the  complete  oxidization  of  the  adatom  and 
the  formation  of  adatom-O-Ce  bonds.  The  present  results  are  sup¬ 
ported  by  previous  experimental  investigations.  For  example,  Fu 
et  al.  [30]  suggested  that  there  are  stable  Au-O-Ce  structures  in 
Au/Ce02.  Using  the  X-ray  diffraction  (XRD)  and  X-ray  absorption 
fine-structure  (XAFS)  techniques,  similar  argument  can  be  made  for 
the  Cu-Ce02,  Pt-Ce02  and  Pd-Ce02i  in  which  the  M-O-Ce  bonds 
are  formed  [14,15,72]. 

We  find  that  the  order  on  the  adsorption  energy  of  metal  (Cu, 
Ag,  Au)  on  Ce02(l  1 1)  follows  the  Cu>Ag>Au.  This  is  consistent 
with  the  conclusion  in  Ref.  [21].  Our  further  analysis  of  electronic 
structures  indicates  that  the  excess  electron  upon  the  M  (M  =  Cu, 
Ag,  Au)  adsorption  can  indeed  localize  on  the  different  Ce  ions.  For 
the  same  adsorption  site,  a  general  rule  is  that  the  electron  prefers 
to  localize  on  Ce  ion  which  is  far  from  the  adsorbed  M  (M  =  Cu,  Ag, 
Au).  However,  there  are  some  differences  among  Cu,  Ag,  and  Au 
adsorption.  For  the  case  of  Ag  adsorbed,  the  Ou-Ou  and  Ce-Ce  site 
can  spontaneously  converge  to  Ou  site. 

It  must  be  emphasized  that  most  of  configurations  with  the 
excess  electrons  localized  at  different  neighboring  Ce  ions  can¬ 
not  be  derived  from  direct  relaxation  from  the  ideal  Ce02(l  1 1) 
surface.  For  example,  for  the  undistorted  structures  of  Cu  and  Ag 
adsorbed  on  Ou  site,  the  excess  electron  delocalizes  at  three  NN 
Ce  ions  and  structures  are  not  stable  with  the  smaller  adsorption 
energies  by  0.2-0.35  eV  as  compared  with  the  distorted  cases.  Thus, 
structural  distortion  is  necessary  to  achieve  a  desired  as  well  as 
localized  configuration.  In  other  words,  the  local  structural  distor¬ 
tion  enables  the  excess  electron  to  be  localized  at  different  sites.  As 
compared  with  the  delocalized  solutions  based  on  the  undistorted 
structures,  the  variations  of  M-0  and  M-Ce  bond  lengths  caused  by 
the  structural  distortions  are  only  in  the  small  range  of  0.01-0.04  A 
and  0.1 1-0.24  A,  suggesting  that  experimental  realization  is  feasi¬ 
ble  by  doping  or  controlling  the  atmosphere.  Recent  investigation 
[73]  regarding  the  effects  for  Ti,  Zr,  and  Hf  doping  on  CO  and  NO 
adsorption  on  Ce02(l  1  0)  surface  has  shown  that  the  dopants  seek 
and  replicate  its  most  favored  coordination  environment,  leading 


to  a  distorted  atomic  structure  of  the  host  oxide.  These  experi¬ 
mental  findings  correlate  well  with  our  work.  The  combinations 
of  experimental  result  with  our  work  deepen  the  understanding 
of  the  catalytic  properties  of  Ce02  surface  supported  metal  clus¬ 
ters  and  could  also  be  expanded  to  other  important  metal/oxide 
systems. 

The  rationalization  why  there  exist  the  various  electron  local¬ 
ization  configurations  can  be  interpreted  by  adopting  the  opinion 
of  Madelung  potential  as  suggested  in  Ref.  [52],  which  states  that 
the  more  energetically  favorable,  the  lower  Madelung  potential. 
Besides,  the  possibility  of  the  direct  f-f  hopping  for  Ce  ions,  the 
geometric  structure,  and  the  intrinsic  characteristics  of  4f  elec¬ 
trons  have  also  been  suggested  to  explain  the  charge  distribution 
in  ceria-based  catalyst  [49,74,75]. 

5.  Conclusion 

Adsorption  mechanisms  of  metal  (Cu,  Ag,  and  Au)  onCe02(l  1 1) 
surface  have  been  investigated  using  the  first-principles  density- 
functional  theory  calculations  considering  the  on-site  Coulomb 
interaction  within  the  PBE  +  U  scheme.  The  systematical  analy¬ 
sis  of  the  adsorptions  energies,  atomic  and  electronic  structures 
reveals  that  the  energetically  favorable  configuration  for  Cu  and 
Ag  adsorption  is  at  the  3-fold  Od  site,  whereas  one  for  Au  adsorp¬ 
tion  is  at  the  Ou-Ou  site.  Cu  atom  exhibits  the  largest  adsorption 
energy,  followed  by  Ag  and  then  Au.  Adsorption  is  accompanied 
by  a  charge  transfer  process,  which  makes  the  adsorbed  metal 
atom  completely  oxidized  and  one  of  the  neighboring  Ce4+  reduced. 
Charge  redistribution  with  electron  localized  on  different  Ce  ions 
neighboring  to  adsorbed  atom  is  observed.  Adsorption  configura¬ 
tion  with  the  excess  electron  localized  on  the  surface  Ce  ion  away 
from  the  adsorbed  metal  is  energetically  more  favorable.  It  has  also 
been  found  the  structural  distortion  play  a  decisive  role  for  under¬ 
standing  metal  atoms  adsorbed  on  Ce02(l  1 1 )  surface.  The  present 
methodology  and  observed  phenomena  are  not  only  significant  for 
theoretical  investigations,  but  also  provide  a  rational  guidance  for 
the  experimental  studies. 
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